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Abstract. Adopting neutrino oscillation parameters obtained by Super-Kamiokande, 
a numerical computer experiment for neutrino events occurring outside the detector, 
is carried out in the same SK live days, 1645.9 live days, constructing the virtual 
Super-Kamiokande detector in the computer. The numerical results by the computer 
experiment could be directly compared with the real SK experimental data. The 
comparison shows that it is difficult to obtain convincing conclusion on the existence 
of the neutrino oscillation with specified neutrino oscillation parameters claimed by 
SK through analysis for neutrino events occurring outside the detector. 
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1. Introduction 

The present situation around — v T oscillation initiated by Super-Kamiokande - 
hereafter, SK, simply — is to be described as follows if we follow majority opinion: 

(1) SK found anomaly in the atmospheric neutrinos (muon deficit) [T], 

(2) SK found anomaly in the zenith angle distribution [2], 

(3) SK discovered atmospheric neutrino oscillation [3], 

(4) Soudan 2 and MACRO confirmed SK results @], 

(5) K2K confirmed atmospheric neutrino results based on the baseline experiment j5], 

(6) SK found oscillations themselves directly [5J, 
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(7) MINOS begins precise experiment and confirmed SK results [7]. 

As one fully understand from the history on — v T oscillation on atmospheric 
neutrinos, it is said that the experimental results by SK play decisively important role 
in — v T oscillation. Therefore, to reach the final conclusions on — v T oscillation by 
SK, it is more desirable that careful examination on the SK conclusion is performed by 
the minority opinion. 

We have been performing Computer Numerical Experiments to re-examine the 
atmospheric neutrino experiment results, constructing the Virtual Super-Kamiokande 
apparatus in the computer and analyzing the virtual neutrino events produced both 
inside and outside the detector. In principle, we could obtain any physical events which 
SK obtain and, therefore, could compare our results with SK results directly, which may 
lead less ambiguous discussion around the atmospheric neutrino oscillation problems. 

In the discussion on the atmospheric neutrino oscillation, there are two fundamental 
issues to be clarified from the experimental point of view, namely, the discrimination 
of the neutrino events (electron events or muon events) and the differences between the 
directions of the incident neutrinos and those of the produced leptons. 

SK treats two different types of the neutrino events as for the points of their 
occurrence with regard to the detector. Ones are the neutrino events occurring inside 
the detector, namely, Fully Contained Events and Partially Contained Events in the SK 
terminology and the others are neutrino events occurring outside the detector, namely, 
Upward Through Going Muon Events and Stopping Muon Events in their terminology. 

The mutual relation between two fundamental issues mentioned above and different 
types of the neutrino events is as follows: in the neutrino events occurring inside the 
detector, one should be careful for the discrimination between muon neutrino events and 
electron neutrino events, because muon neutrino events produce the different Cherenkov 
light patterns from those of electron neutrino events inside the detector. Furthermore, 
one should be very sensitive to the differences between the directions of the incident 
neutrinos and those of the emitted leptons ( muons or electrons) , because their scattering 
angles could not be neglected due to their smaller energies. 

On the other hand, in the neutrino events occurring outside the detector, one need 
not worry about the both discrimination among the neutrino events and scattering 
angles of the neutrino events, because these events should be identified as the muon 
events due to their larger effective volume and the scattering angle of the produced 
leptons concerned could be neglected due to their higher energies. 

In the present paper, we analyze the neutrino events occurring outside the detector. 
The reason why we carry out the analysis of Upward Through Going Muon Events and 
Stopping Muon Events is that this analysis is more simpler and easier than that of the 
Fully Contained Events and Partially Contained Events. Our analysis of the neutrino 
events occurring inside the detector will be presented as a series of the papers elsewhere. 

| SK claim they discriminate muon neutrino events from electron neutrino events almost perfectly [8J, 
while Anokhina and Galkin [9] and Galkin et al[10] have different opinions. 
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2. The Algorithm of the Computer Numerical Experiment 

Our fundamental approach to examine the experimental results obtained by SK is to 
carry out "re-experiment" of SK. For the purpose, we have constructed the virtual SK 
detector in the computer and analyze virtual "experimental" events produced in the 
virtual detector. 

The characteristics of our Computer Numerical Experiment are that we could 
obtain any physical results which SK obtain and we could examine SK results more 
concretely, comparing SK results with our results directly. 

Here, we give the algorithm for our computer numerical experiment. In Figured! we 
show the concept of our numerical computer experiment schematically. The procedure 
of our numerical experiment is as follows: we start a neutrino with some energy on 
the opposite side to the Earth to the Super-kamiokande. For the purpose, taking the 
atmospheric neutrino energy spectrum on the opposite side of the Earth to the detector 
which SK utilizes [TT] , we choose a neutrino with some energy by the random sampling 
from this spectrum. 




Distance from the detector (km) 

Figure 1. Schematic illustration Figure 2. Densities of the Earth 

of the numerical experiment. for different zenith angles 6 V . 

Neutrinos thus chosen penetrate into the interior of the Earth and interact at some 
point inside the Earth after the traverse through the layer with different densities 
shown in Figure [2] [12] . The interaction points of the neutrino events concerned are 
determined by the random sampling, following the probability function for the given 
neutrino interaction mean free path where the change of the density of the Earth is taken 
into account, referring to the Preliminary Reference Earth model [12]. The interaction 
mean free paths of the neutrinos for given energies concerned depend on the both density 
of the location and the zenith angle of the incident neutrino. 

In Figure [3] and Figure HI we give them in cos 6^ = —1.0 and cos6 v = —0.7, 
where cos 9 V denotes zenith angle of the incident neutrino, as examples. Furthermore, 
the neutrino concerned produces a daughter muon by deep inelastic scattering [3], the 
energy of which is also determined by the random sampling, following the probability 
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function for the charged current (deep inelastic scattering) differential cross section 
concerned [T3] . 

As for Upward Going Muon Events and Upward Stopping Muon Events in our 
computer numerical experiment, we could get the essential information around them, 
namely, what nature the neutrinos concerned have ( normal or anti-), their location 
where they are produced and what energies they have, what category they belong to 
( Upward Through Going Muon Events or Upward Stopping Muon Events ), while the 
real experiments by SK give only the discrimination between Upward Through Going 
Muon Events and Upward Stopping Muon Events. 
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Figure 3. Mean free paths of the 
neutrino as the function of energy 
for zenith angle (cos# = —1.0) for 
different energies. 
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Figure 4. Mean free paths of the 
neutrino as the function of energy 
for zenith angle(cos# = —0.7) for 
different energies. 



Exactly speaking, the procedure which we really utilize is slightly different from 
what we just mention above, although it is mathematically equivalent to what we 
mention above. 

Our real procedure adopted is as follows: instead of sampling the neutrino concerned 
directly from the incident neutrino energy spectrum on the surface of the Earth and 
pursing it from its interaction point in stochastic manner, we define the neutrino 
interaction probability function at t n which is measured from the surface of the Earth, 
in the following. 

, , , dt \ / dt 

Pint{K(T>),tn, COS 9 v )dt 



n i 



1=1 



\{E v (p), ti, p{9 v , U)) J \\ n (E„(p),t n , p{Q v ,ti)) 



1) 



, where \%(E v (m, ti, p{0 v , ti)) is the mean free path of the neutrino interaction for normal 
and anti-neutrino with E v (v) at the depth ti, respectively where the density of the Earth 
is p{Q v ,ti) and cos^ denotes nadir angle. Here, we consider the deep inelastic charged 
current interaction [13] and density of the Earth is taken from the Preliminary Reference 
Earth Model [12] . Here, we adopt the deep inelastic scattering exclusively as the cause 
of neutrino interaction, because other neutrino interactions could be neglected for our 
case. 
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By combining Eq.flTJ with the atmospheric neutrino energy spectrum at the 
opposite side of the Earth to SK, we construct N intjnu u(E u (p), t iy cos 9 v )dl, neutrino 
interaction energy spectrum without neutrino oscillations for normal and anti-neutrino 
which denotes frequency of the neutrino interaction in (t n , dl) in the following, 

N int ,nuu(Eu(v), t n , cos9 u )dl = N sp (E v{v) , cos9 u )dl x P int (E v ( V) ,t n , cos9 u ) (2) 

, where N ap (E v ^, cos9 v ) denotes the atmospheric neutrino spectrum for normal and 
anti- neutrino at the opposite side of the Earth to SK detector 

The neutrino interaction energy spectrum, N inttnu u(E v (p^ , t n , cos 9 v )dl, , denotes 
frequency for neutrino interaction in (t n , dl)for null oscillation. The dl chosen in relation 
to the maximum energy of neutrino energy spectrum in the following. 

The maximum energies of the emitted muons are equal to energies of the parent 
neutrino with maximum primary energy. We utilize the incident neutrino energy 
spectrum by Honda et al. [11] which is also utilized by SK, the maximum energy of 
which is 10 TeV. Therefore, we adopt 10 TeV as the maximum energy of muon to be 
pursued in our computer numerical experiment. 
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Figure 5. Range Straggling of Figure 6. Range Straggling of 1 

10 TeV muon in the standard TeV muon in the standard rock, 

rock. 

Now, we decide the maximum range for 10 TeV muon which is necessary for our 
computer numerical experiment. For the purpose, we simulate muons with 10 
TeV in stochastic manner exactly, taking into consideration of the all processes of 
bremsstrahlung, direct electron pair production, nuclear interaction and ionization loss. 

In Figure we give the range fluctuation of 10 TeV muon thus obtained in the 
standard rock. It is easily understood from the figure that the range of the muon 
is widely distributed by the fluctuation due to nature in these elementary processes. 
The shorter ranges than their average are caused by catastrophic energy loss due to 
the bremsstrahlung and nuclear interaction (see, for example, Figure [TT1) . while longer 
ranges than their average are caused essentially by lager number of the direct electron 
pair production ( continuous energy loss-like) (see, for example Figure fT2j) . From Figure 
[5j we could adopt 2.0 x 10 6 g/cm 2 as the maximum range of 10 TeV muon. In other 
word, when muons with 10 TeV are produced at the place far 2.0 x 10 6 g/cm 2 from 
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the detector, such muons never reach the detector. Therefore, it is enough for us to 
have interest in high energy neutrinos only which are produced within 2.0 x 10 6 g/cm 2 
from the detector. Now, we take 2.0 x 10 6 g/cm 2 as dl which is defined in Eq.(J2]). The 
generation points of such muons distribute uniformly within dl, because dl is smaller 
by higher order magnitude of the mean free path of the neutrino for deep inelastic 
scattering. In Figure [HJ we give the range straggling with 1 TeV muon obtained in 
similar manner as in Figure [51 From Figures [5] and [6], it is easily supposed that the 
effect of the range fluctuation in higher energy muon could not be neglected in spite of 
the rapid decrease of the higher energy muon in muon energy spectrurrj§| . 

Our computer numerical experiment for definite zenith angle of the incident 
neutrino is being carried out as follow: 

Procedure A: We calculate the neutrino interaction probability from Eq.([T]). The value 
of t n is defined as the distance from the surface of the Earth to dl. Next, we calculate 
the interaction neutrino energy spectrum from Eq.(j2J). 

Procedure B: By using N int (E u (p), t n , cos 6 u (p))dl thus obtained from Eq.©, we sample 
randomly E v (v), the energy of the incident (anti-) neutrino to be pursued in the following 
way. 

s = / N intinu ii(E u (u), t n , cos 9 v ^f))dE v (p) I I N intjnu u(E v ^,t n , cos 9^)) dE^ 

■^v(y),min j ^u(y) ,min 

, where E v (p^ itnax , E v ^^ min denote the maximum and minimum energy of the neutrino 
concerned, respectively. And £ denote uniform random number between (0,1). 
Procedure C: We determine the interaction point of the (anti-) neutrino event with 
energy determined by Procedure B. By using £, the uniform random number between 
(0,1), we determine the interaction point £dl. Namely, the neutrino interaction occurs 
from the detector, t« is given as, 

U = {l-i)dl (4) 

The detailed study on the pursuit of high energy muon, relating to the range energy 
relation, will appear in a separate paper elsewhere. 

Procedure D: The E^), the energy of the emitted (anti-)muon which occur at (1 —£)dl 
from the detector produced by E v is determined from the following equation, 

r^fiQI) j P ^ jj,CJT) .max 

id= J _ D^E^^E^jdE^ J ^ D(E u(P ), E^jdE^ (5) 

, where D(E v (p), E^-p^dE^ is the differential cross section of the deep inelastic 
scattering for (anti-) neutrino [13] and E^(ji^ max equal to the E u (p\. The £d is a uniform 
random number between (0,1). 

Procedure E: For the (anti-)muon whose energy and production point is determined 
from Procedures C and D, we examine the behavior of the trajectory of the (anti-)muon 
toward the SK detector in a stochastic manner. Namely, each individual muon is pursued 

§ The range energy relation of high energy muon will be studied in detail in a separate paper elsewhere. 
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For given cos V (= cos 0m) 



Random sampling of Ev from the neutrino interaction energy spectrum 
at Kamioka cite( Eq.(3)(without oscillation) and Eq.(8)(with oscillation) 

T 



Random sampling of E^ from the deep inelastic scattering 
differential energy spectrum (Eq.(5)) 



I 



Random sampling of (1- £)dl, the interaction point 

of the neutrino event concerned from the detector (Eq.(5)) 



T 



Random sampling of the interaction point (x,y) in the two dimensional plane 
which has one to one correspondence with regard to the SK detector 
for given cos0 v (= cos 0^ (Eq.(4)) 



T 



The pursuit of the muon concerned whose energy is E and the interaction 
point ( x, y, (1- i; )dt ) and goes foward the detector in the exact stochastic 
manner ( Fig. 10, Fig. 11, Fig. 12 ) 
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Figure 7. Grand structure for the computer numerical experiment 



exactly by taking into consideration of the processes of bremsstrahlung, direct electron 
pair production, nuclear interaction and ionization loss in the stochastic manner without 
utilizing the average behavior of the muon concerned, the procedure of which is adopted 
by SK [14] . As the result of it, we determine finally, which category each individual 
muon pursued falls into: [a] stopping before it reaches the detector, [b] stopping inside 
detector, or [c] passing through the detector. 

We repeat Procedures A to E and obtain the neutrino events without oscillation for a 
given live time for the real experiment. In our computer numerical experiment, we have 
accumulated the events concerned which correspond to the real live time for SK. For each 
neutrino event, we know E u ^, the energy of the parent neutrino, E^), the energy of the 
daughter (anti-)muon,(l — £)dl, the interaction point, cos Q v (j>) ( cos ^(m)), the direction 
of both the incident (anti-)neutrino and the emitted (anti-)muon, respectively, and t n , 
the distance between the interaction point of the neutrino events and the opposite side 
of the Earth. 

In the case without neutrino oscillation, we repeat from Procedures A to E, starting 
from the interaction energy spectrum(Eq.([2])). Finally, we obtain frequencies for Upward 
Stopping Muon Events and Upward Through Going Muon Events for given zenith angle 
for null oscillation. Thus, we obtain the zenith angle distributions for these categorically 
different events without neutrino oscillation after summation over zenith angle. 
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In the case with neutrino oscillation, Ni nt {E u (p) : t n , cos8 v (p-))dt in Eq.([2]) should be 
replaced by 

N in t,osc(E v ^,t n , cos 9 v (p))dt = N inttnuU {E v{V ),t n , cos B v{v) )dt x P(i/ m (f^) — ► ^(%)) (6) 
, where 

r,/ / x /xx 9«« 9 /l.27Am 2 (eK 2 )t„(£;m)\ 

, the survival probability of a given v^. 

Here, we adopt sin 2 26 = 1.00 and Am 2 = 2.1 x 10 _3 eV 2 [TH] , as obtained from 
SK experiment. Other procedures are exactly same as in the case without neutrino 
oscillation. 

Thus, Eq.(j3]) should be replaced by 

Ni n t,osc{E u (p), t n , cos v (p})dE v fp) I I N inti0SC (E u (y-),t n , cos 9 v ^)dE v ^ 

-^vCv) ,min / .min 

, by using Eq.(J5]) and other procedure are not changed. 

The Procedures A to E in both cases without oscillation and with oscillation is 
illustrated in Figure [71 

3. Results and discussion 

3.1. The differential neutrino interaction energy spectrum with neutrino oscillation 

In Figure [U we give the differential neutrino interaction energy spectrum with oscillation 
(Eq.(jSD) for different zenith angle from which we randomly sample the energy of the 
neutrino event concerned. 

From the figure, it is easily understood that the differential neutrino interaction 
energy spectrum are not so influenced by the characteristics of the nature of oscillatory 
above ~ 40 GeV under SK neutrino oscillation parameters. The result of the random 
sampling by Eq.flB]) is given in Figure EB The excellent agreement between the theoretical 
curve (solid line) and the results sampled (histogram) guarantees the validity of the 
sampling procedure by Eq.(jHJ). 

The nature of oscillatory in the neutrino interaction energy spectrum a play crucial 
role below ~ 40 GeV under the SK neutrino oscillation parameters, on which we will 
examine in a series of the subsequent papers in detail. 

3.2. Discrimination between Upward Through Going Muon Events and Stopping Muon 
Events 

The physical origin of Upward Through Going Muon Events is the same that of Stopping 
Muon Events and the only difference between them lies in the situation that the former 
may eventually pass through the detector, while the latter may eventually stop the both 
initiated by the deep inelastic scattering. 

The generation point of the neutrino event concerned is given Eq.(T4]). As the 
generation point of the neutrino event distribute over dl (2.0 x 10 16 g/cm 2 in our case 
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) randomly shown in Figure dOj the most biggest fluctuation in the muon behavior 
essentially comes from the generation point of the neutrino events over dl. Such 
fluctuation could not be considered in the Detector Simulation by SK. 

The energy of the neutrino is sampled randomly from Eq. (jSJ) and the energy of the 
emitted muon from the neutrino event concerned is sampled randomly from Eq. (|5]) . The 
muons whose energies and whose generation points thus determined are simulated in the 
stochastic manner exactly, taking into account of the all processes of the direct electron 
pair production, bremsstrahlung, nuclear interaction and ionization loss. Namely, these 
muons thus sampled lose their energies in stochastic manner in the course of their 
passages toward SK detector. 

In summarized speaking, the fluctuation effect in the muon behavior secondarily 
comes from the range fluctuation of the muon, although this effect is rather small 
compared with the effect from the generation point of the neutrino events but could 
not be neglected. 

Even if the physical natures of the muons are exactly same, some muon may be 
classified as the Upward Through Going Muon Events or may be done as Stopping 
Muon Events eventually due to the complicated geometry — configuration of (x, y) - 
for muons concerned, which is schematically shown in Figure [101 

We take into account of the complicated circumstances for the muons concerned as 
exactly as possible. As the results of it, we could discriminate finally Upward Through 
Going Muon Events from the Stopping Muon Events. 

Among elementary processes contributed to muon energy losses, the process of the 
direct electron pair production, bremsstrahlung and nuclear interaction play important 
roles from the point of fluctuation. In Figures [11] and [12], how these three elementary 
processes take place randomly in the stochastic ways are indicated as the functions of 
the traversed distance in standard rock by a muon with 10 TeV. Therefore, we treat the 
muon energy loss in stochastic manner exactly. 
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Figure 10. The distribution of the generation point of the neutrino event 
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Figure 11. An example of a real 
high energy muon behavior. 



Figure 12. An another example 
of a real high energy muon 
behavior. 



In the figures, we show graphically where, which interaction (electron direct pair 
production, bremsstrahlung or nuclear interaction) occurs and what energy with respect 
to the primary energy is lost by the interaction concerned. The /x denotes the sustained 
energy of the muon concerned as the result of all possible energy losses. The D denote 
the energy loss due to direct electron pair production, B is the energy loss due to the 
bremsstrahlung and N denotes nuclear interaction, respectively. The Ratio denotes the 
ratio of energy loss concerned to the primary energy. 

In Figure [TTJ we give an example of a muon with shorter range, say, 1500 meter. 
We could understand the muon loses about 30 % catastrophically at the early stage by 
the bremsstrahlung and bigger energy losses may come from the bremsstrahlung and 
nuclear interaction and smaller energy losses essentially come from the direct electron 
pair production. In Figure [12] which is contrast to Figure [TTJ we give an example of a 
muon with range over 3000 meters. In this case, bremsstrahlung and nuclear interaction 
do not contribute to bigger energy losses and the accumulation of large number of the 
direct electron pair production events with rather smaller energies result in bigger energy 
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loss. 

On the contrary to our exact treatment on the individual muon behavior, SK seem 
to treat them in the sense of as the average theory [TH [T5] . SK utilize the following 
conventional formula, revised by Lohmann, Knop and Voss [16J. 
dE 

— = a + bx E (9) 

ax 

Then, this description expresses the average energy loss of the muon which does not 
give the range straggling of the muon and dose not expresses the real energy loss by 
individual muon. When the energy of the muon is not so pretty high where we could 
neglect the fluctuation effect of the muon range in individual case, we could apply this 
expression to the analysis of real behavior of the muon. However, the energy region for 
muon where we are interested in is now rather high and, therefore, we should pursue 
the muon behavior in the stochastic manner exactly, but in not average theory, even if 
the fluctuation effect is taken into account [jj. It should be emphasized that the average 
picture of the muons obtained from the average theory like Eq.Q does not describe the 
individual behavior of the muon concerned and the average picture of the muons is the 
average behavior of the individual muon which is deviated from their average behavior 
when we could not neglect the fluctuation effect coming from their stochastic nature 
in the elementary processes. Furthermore, it should be noticed that the stochastic 
treatment to the individual muon becomes important in the analysis of statistically 
insufficient number of the events concerned. 

In summarized speaking, the fluctuation effect in the muon behavior secondarily 
comes from the range fluctuation of the muon, although this effect is rather small 
compared with the effect from the generation point of the neutrino events but could 
not be neglected. 



|| Recently, SK introduce the concept of showering muon in their treatment on the high energy upward 
muon event [17j . In their paper, they seem to treat high energy muon in stochastic manner, utilizing 
GEANT3. 



The numerical computer experiment for the neutrino events 



edited Feb. 08, '08 12 



3.3. Zenith Angle distribution for Upward Through Going Muon Events 



In Figure [13] we compare the experimental results by SK with corresponding our 
computer numerical results for both with neutrino oscillation and without neutrino 
oscillation in the case of Upward Through Going Muon Events for the same 1645.9 live 
days |15j . From the figure, we could not conclude the experimental data agree with null 
oscillation or with oscillation. 

The 1645.9 live days is not enough for extracting definite conclusion on the neutrino 
oscillation under SK neutrino oscillation even if it is possible to get SK neutrino 
oscillation parameters through the analysis of Upward Through Going Muon Events, 
because the statistical fluctuation from the average values is not small. If we compare 
our results with oscillation to our results without oscillation in the figure, we could find 
the event number with oscillation is larger than that without oscillation in the three bins 
(0.0 ~ -0.1), (-0.1 ~ -0.2) and (-0.9 ~ -1.0) . Such phenomena come from the statistical 
fluctuation due to insufficient sampling (accumulation) of the neutrino events. 
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Figure 13. Zenith angle distri- 
bution for the Upward Through 
Going Muon Events for 1645.9 
live-days. 



Figure 14. Zenith angle distri- 
bution for the Upward Through 
Going Muon Events for 164590 
live-days, which is normalized to 
1645.9 live-days. 



In Figure [TU we give corresponding results from our computer numerical experiment for 
164590 days, one hundred times of the SK real live days. Now, we understand from the 
figure that our results with oscillation are always smaller than that without oscillation 
which attributes to larger statistics. This denotes that sampling numbers for 164590 live 
days is quite enough for extracting the definite conclusion on the neutrino oscillation 
around the SK neutrino oscillation parameters. 

Furthermore, it is reasonable to say from the figure that we could not discriminate 
the difference in both with null oscillation and with oscillation under SK oscillation 
parameters from our computer numerical experiment. In other words, this is quite 
understandable, because the definition of the Upward through Going Muon Events is the 
muon events, the minimum energy of which is only given ~ 1.6 GeV and nothing more 
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both the incident energy, say, nobody knows both generation points and termination 
points and energies of the events concerned, but knows their zenith angle merely. 
It should be rather curious, if we really could extract definite neutrino oscillation 
parameters from such extremely poor information around the events concerned. 

In Figure HU we add, on purpose, the SK experimental data for 1645.9 live days 
together with our data for 164590 live days — 450 live years — as one hundred times as 
SK real live days. From the comparison between SK results for 1645.9 live days and 
our result for 164590 live days — 450 live years — , we could not conclude directly that 
SK experimental data rather agree with our results with oscillation, because statistics 
of both results is quite different from each other and we do not allow to compare them 
directl}0 . 

Comparing the smoothness of the histograms in Figure [13] with these in Figure 
[T4"l we could conclude that the average values of zenith angle distributions attain 
sufficiently at the 'true' ones. In other words, we think SK live days does not give enough 
statistics to extract clear cut conclusion under SK neutrino oscillation parameters from 
the analysis of Upward Through Going Muon Events, even if SK finally could give 
definite neutrino oscillation parameters from the analysis of Fully Contained Events 
and Partially Contained Events. 

In Figure [T5l we give the correlation diagram between the energies of the emitted 
muons and the distance from their occurrence points up to the detector for Upward 
Through Going Muon Events for null oscillation, pursuing the muons concerned in 
stochastic manner. The average energy of the emitted muons is 390 GeV and their 
average distance 161 meter. 

In Figure [T6J we give the corresponding quantities for oscillation. The average 
energy of the emitted muons is 402 GeV and their average distance is 168 meters. It is 
easily understood from the Figure [15] and Figure [16] that the contribution from higher 
energy part of the muon could not negligible, which is larger beyond the expectation by 
SK whose analysis is based on the average theory for muon propagation. Furthermore, 
comparing the Figure [15] with Figure [16] we could not find clear differences between 
them, which denote the difficulty in extracting the SK oscillation parameters from the 
analysis of the Upward through Going Muon events. 

% SK discusses the evidence of the neutrino oscillation under SK oscillation parameters, based on the 
comparison of the real experimental data (for example, for 1645 live days) with Monte Carlo Simulation 
data without oscillation (for example, 100 live years which are much larger than that of real experiment). 
However, it seems not suitable to compare directly the one with smaller statistics with the other with 
further large statistics, because we have the possibility to kill real fluctuation. SK is desirable to 
compare the real experimental data for 1645.9 live days with Monte Carlo simulation for the same live 
days to extract definite conclusion. Further more, it should be noticed that a 1645.9 live days might 
be different from another 1645.9 days due to significant insufficient statistics. 
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Figure 15. 



Figure 16. 



Figure 17. 



Figure 1151 Correlation diagram between the emitted muon energies and their 
distance from their generation points to the detector for null oscillation in the case 
of Upward Through Going Muon Events. 

Figure 1161 Correlation diagram between the emitted muon energies and their 
distance from their generation points to the detector for the oscillation in the case 
of Upward Through Going Muon Events. 

Figure II TL The distance distribution for Upward Through Going Muon Events in 
both cases of null oscillation and the oscillation. 

In Figure [T7J we give the distance distribution for Upward Through Going Muon Events 
in the cases of both null oscillation and the oscillation for 1645.9 SK live days. The 
distance denotes the interval from the generation point of the events concerned to the 
detector. It is easily understood from the figure that we could not discriminate the 
difference in both cases well, which suggests difficulty in extracting the definite SK 
neutrino oscillation parameters. 
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Figure 18. Comparison of 
the zenith angle distribution for 
Upward Through Going Muon 
Event by SK with their Monte 
Carlo Results 



Now, we return the SK original data for Upward Through Going Muon Events which 
we reproduce it in Figure [THJ Now, we point out clear contradiction in Figure [18j The 
first to be examined is that frequency of the experimental data is significantly larger 
than that of their Monte Carlo Simulation without oscillation in the wider region (-0.6, 
0.0) for cos6 and the second is rather nice agreement between them in the region (- 
1.0,-0.6) for cos 6 . The first denotes that usually we should expect smaller frequency 
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with oscillation than that without oscillation. However, here, the situation becomes 
reverse which contradict with the commonsense. Furthermore, in the region (1.0,-0.6) 
for cos#, the effect of the neutrino oscillation should be pronounced due to longer 
path length of the incident neutrino. However, the real situation seems to indicate no 
oscillation. In this case, there are three possible interpretations to such a situation. 
Namely, [a] Monte Carlo Simulation is not unreliable, while the experimental data is 
reliable including statistical uncertainty, [b] Monte Carlo Simulation is reliable, while 
statistics of the experimental data is too small to give larger frequency of the neutrino 
events compared with Monte Carlo events, [c] Neither experimental data nor Monte 
Carlo data is unreliable to exclude the definite oscillation parameters. 

Furthermore, from the view point of the comparison between experimental data 
and Monte Carlo one, we should make two comments: the fist is that SK Monte Carlo 
simulation does not consider the real behaviors of the muons in stochastic manner, but 
their average behaviors only. The second is that it is not suitable way to compare 
experimental data for 1645.9 days with Monte Carlo one for 100 live years, because of 
big difference in statistics between them. The normalization of 100 live yeas to 1645.9 
days never guarantees the validity of their transformation without examination on the 
validity concerned. 

To clarify the contradiction in Figure [THJ, it is necessary for SK to disclose their 
Monte Carlo procedure in detail. 



3.4- Zenith angle distribution for the Stopping Muon Events 
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Figure 19. Zenith angle dis- 
tribution for Upward Stopping 
Muon Events for 1645.9 live-days 



Figure 20. Zenith angle dis- 
tribution for Upward Stopping 
Muon Events for 164590 live- 
days, which is normalized to 
1645.9 live-days. 



It is natural to think that Upward Stopping Muon Events are much influenced by 
fluctuation, compared with Upward Through Going Muon Events, because they are 
more directly influenced by the complicated geometry of the detector with regard to the 
incident neutrinos than Upward Through Going Muon Events are done. 
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In Figure dHJ we compare the SK experimental data for 1645.9 live days with our 
data with oscillation and without oscillation for Upward Stopping Muon Events for the 
same 1645.9 live days. From the figure, also, we could conclude that SK experimental 
data neither agree with our data with oscillation nor that without oscillation. In the 
figure, it should be noticed that the number of the neutrino events with oscillation 
is larger than that without oscillation due to fluctuation effect for cos6> = —0.05 
(0.0 0.10) . 

In Figure [201 we gi ye our results with oscillation and without oscillation for 164590 
live days together with SK experimental data which, in principle, should not be allowed 
to be compared with our data due to big difference in statistics as stated in the case of 
Upward Through Going Muon Events (see, also footnote on page 13). Compared Figure 
20l with Figure [19] we understand that the difference between quantity with oscillation 
and that without oscillation in Figure [20] decrease compared with that in Figure [19] and 
both histograms in Figure [20] become smooth compared with those in Figure [19] which 
is evidence that "average value " attain at ideal average value. 
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Figure 21. Figure 22. Figure 23. 

Figure 1211 Correlation diagram between the emitted muon energies and their 
distance from their generation points to the detector for null oscillation in the case 
of Upward Stopping Muon Events. 

Figure 1221 Correlation diagram between the emitted muon energies and their 
distance from their generation points to the detector for the oscillation in the case 
of Upward Stopping Muon Events. 

Figure 1231 The distance distribution for Upward Stopping Muon Events in both 
cases of null oscillation and the oscillation. 

In Figure we give the correlation diagram between the emitted muon energies and 
their distances in Stopping Muon Events for null oscillation. The average energy of the 
emitted muons concerned is 55 GeV and their average distance is 44 meters. In Figure 
22l we give the similar diagram for SK oscillation parameters as Figure [21] The average 
energy of the muons is 56 GeV and their average distance is 45 meters. The situation 
around Figure [21] and Figure [22] is essentially same as in that around Figure [15] and [16] 
In Figure [23] we give the distance distribution for Upward Stopping Muon Events 
in both cases of null oscillation and the oscillation. It is reasonable to say from the 
figure that we could not find the significant difference in the both distance distribution, 
which may lead to the difficulty in extracting the definite oscillation parameters like SK 
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do. 



3. 5. The energy spectrum for the parent neutrinos and their daughter muons in the 
presence of the neutrino oscillation 
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Figure 24. Figure 25. Figure 26. 

Figure I241 Energy spectra for parent neutrinos and daughter muons for Upward 

Through Going Muon Events for cos 9 V = —0.957 (—0.95 ~ —1.0). 

Figure I25L Energy spectra for parent neutrinos and daughter muons for Upward 

Through Going Muon Events for cos6 u = -0.525 (-0.500 ~ -0.550). 

Figure 1261 Energy spectra for parent neutrinos and daughter muons for Upward 

Through Going Muon Events for cos 9 U = —0.025 (—0.0 ~ —0.05). 
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Figure 27. Figure 28. Figure 29. 

Figure I27L Energy spectra for parent neutrinos and daughter muons for Upward 
Stopping Muon Events for cos6 u = —0.957 (—0.95 ~ —1.0). 

Figure I28L Energy spectra for parent neutrinos and daughter muons for Upward 
Stopping Muon Events for cos9 u = —0.525 (—0.500 ~ —0.550). 
Figure 1291 Energy spectra for parent neutrinos and daughter muons for Upward 
Stopping Muon Events for cos^ = —0.025 (—0.0 ~ —0.05). 

Now, we examine energy spectra of the parent neutrinos and their daughter muons for 
the Upward Through Going Muon Events and the Stopping Muon Events with different 
cos 6^(6^). 

In Figures EH [25] and [26J, we give energy spectra of Upward Through Going Muon 
Events for cos^ = —0.975, (upward) — 0.525(diagonally) and — 0.025(horizontally), 
respectively for SK live days, 1645.9 days. In Figure [271 [28] and [29], we gi ye the energy 
spectra of the Stopping Muon Events for the same for the same SK live days. It seems 
to be clear from Figures EH to [29] that all these energy spectra are never smooth which 
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denote insufficient statistics to extract definite conclusions, particularly as shown in 
Figure [26] 

In Figures EH [25] and [26] we give the average energies for parent neutrinos and 
daughter muons ( in parentheses) in the case of the Upward Through Going Muon 
Events for cos0„ = -0.975, -0.525 and -0.025, respectively, namely 529GeV (343GeV), 
657GeV (500GeV) and 658 GeV (407GeV). The corresponding average energy for the 
parent neutrino in SK is ~100 GeV [H] which is pretty smaller than that our values. In 
Figure [271 EE] and [29l we give the average energies of the Stopping Muon Events for the 
same cos^, 290GeV (216GeV), 134GeV (54GeV) and 28GeV (16GeV), respectively. 
The corresponding average energy for neutrino by SK is ~ 5 GeV [H] which is also too 
smaller than our corresponding values. 

As for obtaining the average energy of the Upward Through Going Muon Events, 
SK utilize the average theory which is unsuitable, as already mentioned, for the analysis 
of real behavior of high energy muon and, therefore, they could not consider the 
contribution from higher energy muons correctly. 

As for the estimation of the average energy of Stopping Muon events, it seems to 
be curious that their average energy by SK is nearly same as or a little higher than 
that of Multi-GeV Events. Because Multi-GeV Events comes from essentially single 
meson scattering, while Stopping Muon Events come from deep inelastic scattering, and 
the effective energy for single meson scattering is too lower than that of deep inelastic 
scattering (see, Figure [TTl in the SK paper [15]). The difference in the effective energies 
between Upward through Going Muon Events and Stopping Muon Events comes from the 
difference in the effective ranges from the same elementary process, say, deep inelastic 
scattering. Also, SK estimate that generation points up to the detector distributed over 
meters to several hundred meters [14] . Too small estimation on average energy as 
well as the generation points for Upward Stopping Muon by SK might come from the 
utilization of the average theory on energy loss of the muon. 

Here, let us remark the relation between oscillatory nature of the primary neutrino 
interaction energy spectrum and their influence over Upward Through Going Muon 
Events and Stopping Muon Events. It s easily understood from Figures[2l]to[29l referring 
to Figure El that neutrinos contributed to both Upward Through Going Muon Events 
and Stopping Muon Events for all their zenith angles come from pretty high energies 
which belong non-oscillatory part of the incident neutrino interaction energy spectrum 
under the specified neutrino oscillation parameters by SK. 

In other word, we should say that the analysis of Upward Through Going Muon 
Events and Stopping Muon Events are not suitable means for proving the existence of 
neutrino oscillation under SK neutrino oscillation parameters. 
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3. 6. Comparison of Monte Carlo Simulation in enough statistics by SK with the 
corresponding one by our computer Numerical Experiment 

Here, we compare SK Monte Carlo Simulation in the case of null oscillation for enough 
statistics with our corresponding Computer Numerical Experiment. 

Now, we compare SK Monte Carlo results for 100 live years with ours for 450 
live years for both Upward Through Going Muon Events and Stopping Muon Events 
in the case of null oscillation. Strictly speaking, SK results for 100 live years should 
be compared with ours for same live year. Or our results for 450 live years should be 
compared with the SK results for same live years. However, as SK data for 100 live years 
are supposed to attain at the enough statistics which is free from fluctuation, compared 
with that for 1645.9 days, we might be allowed to compare directly SK data for 100 live 
years with our data for 450 live years. 
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Figure 30. Our Monte Carlo 
result and SK Monte Carlo ones 
for null oscillation for Upward 
Through Going Muon Events. 
Both Monte Carlo results are 
normalized to 1645.9 live-days. 



Figure 31. Our Monte Carlo re- 
sults and SK Monte Carlo ones 
for null oscillation for Upward 
Stopping Muon Events. Both 
Monte Carlo results are normal- 
ized to 1645.9 live-days. 

From Figure [30] we understand that general tendency is similar to each other and 
our results give larger value than SK results. It is reasonable, if we take into account of 
the fact that SK neglect the range fluctuation in high energy muon. 

In Figure [311 we give the direct comparison between SK results and our results for 
Upward Stopping Muon Events. Here, SK results give larger values than our results, 
on the contrary to Figure [301 This reversed situation between Figure [30] and Figure 
l3~T1 seems to be curious, if we take into account of the range fluctuation of high energy 
muon. Also, tendencies in the SK histograms in Figure [31] are somewhat unnatural in 
lack in smoothness in spite of huge statistics. 

It is desirable for SK to disclose their essential information on their Monte Carlo 
method for the case of the null oscillation in order to clarify the difference between 
SK's and ours. The techniques for our Monte Carlo method fully disclose in the present 
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paper. As for our Monte Carl Simulation for the neutrino events inside the detector, we 
disclose our procedure in detail in the subsequent papers. 

4. Conclusion 

Thus, the followings should be concluded: 

(a) The fluctuation effect in the high energy muon behavior from the neutrino event 
(due to deep inelastic scattering) primarily comes from the uniformly randomness of 
the generation point of the neutrino events concerned. This effect could not be treated 
in the correct manner in the Detector Simulation by SK. 

(b) Also, the consideration of the range fluctuation in high energy muon is inevitable, 
although this effect is pretty smaller than that from the generation point of the neutrino 
events. We consider this effect exactly in the stochastic manner. On the other hand, 
SK analyzes the behavior of muons in the average theory, namely, in the concept of 
the average energy loss and average range, not treating them in stochastic manner. 
Combined of [a] with [b], we claim that experimental results for the Upward Through 
Going Muon Events and the Stopping Muon Events do not show the evidence under SK 
neutrino oscillation parameters. 

(c) Generally speaking, it is difficult to imagine rationally to extract so definite 
conclusion around neutrino oscillation parameters, say, sin 2 26 and Am 2 through the 
analysis of the Upward through Going Muon Events and Stopping Muon Events. 
Because, there are too few information around such categorized events to extract the 
definite conclusion on the neutrino oscillation parameters. Such categorized events 
provide the direction of the incident neutrino only. Upward Through Going Muon Events 
never provide us the energies of the incident neutrinos, their interaction point and the 
energies of the emitted muons, their stopping points, giving the information on their 
penetration through the detector only. Stopping Muon Events gives the information 
that they stop inside the detector and other information not obtained by them are same 
as the Upward Through Going Muon Events. 

(d) Even if the SK neutrino oscillation parameters exist really, it is difficult to estimate 
them from the analysis of Upward through Going Muon Events and Stopping Muon 
Events. Because, the effective energies to detect the SK neutrino oscillation parameters, 
say, smaller than 40 GeV, do not lie in the energy region which contribute to the 
production of Upward Through Going Muon Events and the Stopping Muon Events. 

(e) By the definition of the Upward Through Going Muon Events and the Stopping Muon 
Events in the SK, the origins of both events are same(deep inelastic scattering) and 
difference between them is that the former average range is longer than that of the 
latter so that the former average energy is higher than that of the latter. However, the 
latter average is so low which is rather near to the average energy due to multi-GeV 
energy events. This is difficult to be understood. Also, the zenith angle distribution for 
Upward through Going Muon events is difficult to be understood in comparison with 
their Monte Carlo simulation with 100 live years. 
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(f) Even if we could obtain the SK neutrino oscillation parameters from the analysis 
of both Upward through Going Muon Events and the Stopping Muon Events, present 
statistics for experimental data obtained by 1645.9 live days is not enough for getting 
the SK neutrino oscillation parameters. 

(g) As the effective energy to detect oscillatory signature of the neutrino oscillation 
parameters lies below 40 GeV under the assumption of the SK neutrino oscillation 
parameters, first of all, one should concentrate to analyze the single ring muon event 
among the Fully Contained Events, the most clear cut event from the quasi elastic 
scattering in GeV energy region. 

Generally speaking, it is rather difficult to extract so definite physical quantities, 
such as, neutrino oscillation parameters from the atmospheric neutrino beams which 
include much uncertainty factors, such as, the incident neutrino energy spectrum. The 
throughout and concentrated analysis of the single ring lepton events are inevitable, if 
we really extract definite parameters from the cosmic ray experiment in spite of the 
existence of much uncertainty factors involved. 

We are now going to complete our works along this line and will present our results 
from the computer numerical experiment, focusing on the analysis of single-ring muon 
events as a series of subsequent papers elsewhere. 
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